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M. 

The large-scale structure in the distribution of galaxies is thought 
to arise from the gravitational instability of small fluctuations in the 
initial density field of the universe. A key test of this hypothesis 
is that superclusters of galaxies in the process of formation should 
generate systematic infall of other galaxies. This would be evident 
in the pattern of recessional velocities, causing an anisotropy in the 
inferred spatial clustering of galaxies. Here we report a precise 
measurement of this clustering, using the redshifts of more than 
141,000 galaxies from the two-degree-field galaxy redshift survey. 
We determine the parameter /3 = 0°^/b = 0.43 ±0.07, where is the 
total mass-density parameter and b is a measure of the 'bias' of the 
luminous galaxies in the survey. Combined with the anisotropy of 
the cosmic microwave background, our results favour a low-density 
universe with SI ^ 0.3. 
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Hubble showed in 1934 that the pattern of galaxies on the sky is non-random -"^ , and 
successive years have seen ever more ambitious attempts to map the distribution of 
visible matter on cosmological scales. In order to obtain a three-dimensional picture, 
redshift surveys use Hubble's law, v — Hor, to infer approximate radial distances to a 
set of galaxies. The first major surveys of this sort took place in the early 1980s^'^''*'^, 
and were limited to a few thousand redshifts, owing to the limited speed of single-object 
spectroscopy. In the 1990s, redshift surveys were extended to much larger volumes by 
a 'sparse sampling' strategy^. These studies''"'^ established that the universe was close 
to uniform on scales above about 100 Mpc {h = i7o/100kms~^Mpc~^), but with a 
complex nonlinear supercluster network of walls, filaments and voids on smaller scales. 

The origin of this large-scale structure is one of the key issues in cosmology. A 
plausible assumption is that structure grows via gravitational collapse of density fluc- 
tuations that are small at early times - but it is vital to test this idea. One important 
signature of gravitational instability is that collapsing structures should generate 'pecu- 
liar' velocities, 5v, which distort the uniform Hubble expansion. We measure a redshift, 
z, which combines Hubble's law with the radial component of these peculiar velocities: 
cz ~ Hor + S^r ■ f . The apparent density held seen in a redshift survey is thus not a true 
three-dimensional picture, but this can be turned to our advantage. The redshift-space 
distortions have a characteristic form, whose detection can both verify the general idea 
that structure forms by gravitational instability, and also measure the density of the 
universe. The present paper presents measurements of this effect, based on a new large 
redshift survey. 

The 2dF Galaxy Redshift Survey 

New-generation large redshift surveys are made feasible by multiplexed fibre-optic spec- 
troscopy, and the most advanced facility of this sort is the 2-degree Field, mounted 
at the prime focus of the Anglo- Australian Telescope^ which allows 400 spectra to 
be measured simultaneously. The 2dFGRS-'^-'^'-'^'^ was designed to use this instrument 
to measure the redshifts of 250,000 galaxies, to a blue magnitude limit of 6, = 19.45 
(corrected for extinction by dust in the Milky Way). The galaxies were selected from 
an updated version of the APM catalogue-*^^, which is based on scans of photographic 
plates taken with the UK Schmidt Telescope. Survey observations began in 1998, and 
should finish at the end of 2001. At the time of writing, redshifts have been obtained for 
141,402 galaxies. The sky coverage of the 2dFGRS consists of strips in the Northern and 
Southern Galactic Poles (75° x 7.5° in the NGP; 75° x 15° in the SGP), plus a number of 
outlying random fields. Coverage is now sufficiently extensive that near-complete thin 
slices through the galaxy distribution may be constructed, as shown in Fig. 1. This 
image illustrates well the median depth of the survey: approximately z — 0.11. Beyond 
this point, the survey is sensitive only to the more luminous galaxies, and the comov- 
ing density falls rapidly. Nevertheless, the survey volume (including other regions not 
shown) is more than adequate for an accurate determination of the statistical properties 
of galaxy clustering. 

Galaxy correlations in redshift space 

The simplest statistical indicator of peculiar velocities in cosmological structure is the 
two-point correlation function, ^(cr, tt). This measures the excess probability over ran- 
dom of finding a pair of galaxies with a transverse separation a and a linc-of-sight sepa- 
ration TT. In an isotropic universe, this function should be independent of direction, but 
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Figure 1 The distribution of galaxies in part of the 2dFGRS, drawn from a 
total of 141,402 galaxies: the slices are 4° thick, centred at declination —2.5° in 
the NGP (left) and -27.5° in the SGP (right). Not all 2dF fields within the slice 
have been observed at this stage, hence there are weak variations of the density 
of sampling as a function of right ascension. To minimise such features, the slice 
thickness increases to 7.5° between right ascension 13. l'^ and 13. 4^^. This image 
reveals a wealth of detail, including linear supercluster features, often nearly- 
perpendicular to the line of sight. The interesting question to settle statistically 
is whether such transverse features have been enhanced by infall velocities. 

this is not true in redshift space. Transverse separations are true measures of distance, 
but apparent radial separations are distorted by peculiar velocities. This redshift-space 
anisotropy should cause two characteristic effects, operating respectively on small and 
large scales. On small scales, random orbital velocities within galaxy groups cause an 
apparent radial smearing, known as 'fingers of God'. Of greater interest is the large- 
scale effect; if cosmological structure forms via gravitational collapse, there should exist 
coherent infall velocities, and the effect of these is to cause an apparent flattening of 
structures along the line of sight. The general existence of redshift-space distortions 
was recognized in the first redshift surveys^'^'^, but the first comprehensive analysis of 
the phenomenon was performed by Kaiser who showed that they could be used to 
measure the quantity 

where Q. is the cosmological mass density parameter and h is the bias parameter that 
relates the relative density fluctuations of the galaxies and of the total mass: 

^ =h — 

P 
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Figure 2 The redshift-space correlation function for the 2dFGRS, ^(cr, tt), 
plotted as a function of transverse (cr) and radial (tt) pair separation. The func- 
tion was estimated by counting pairs in boxes of side 0.2 hr^ Mpc (assuming an 
= 1 geometry), and then smoothing with a Gaussian of rms width 0.5 h-^ Mpc. 
To illustrate deviations from circular symmetry, the data from the first quadrant 
are repeated with reflection in both axes. This plot clearly displays redshift 
distortions, with 'flngers of God' elongations at small scales and the coherent 
Kaiser flattening at large radii. The overplotted contours show model predic- 
tions with flattening parameter (5 = OP'^ /h = 0.4 and a pairwise dispersion of 
cTp = 400 kms-^ Contours are plotted at ^ = 10,5,2, 1,0.5,0.2,0.1. 

The model predictions assume that the redshift-space power spectrum 
(Ps) may be expressed as a product of the linear Kaiser distortion and a radial 
convolution^^: P,(k) = Pr{k) (1 + (3^^)^ (1 + k^a^^^ /2H^)-\ where ^ = k • f , 
and CTp is the rms pairwise dispersion of the random component of the galaxy ve- 
locity field. This model gives a very accurate fit to exact nonlinear simulations"*^^. 
For the real-space power spectrum, Prik), we take the estimate obtained by de- 
projecting the angular clustering in the APM survey"*"^'"*"^. This agrees very well 
with estimates that can be made directly from the 2dFGRS, as will be discussed 
elsewhere. We use this model only to estimate the scale dependence of the 
quadrupole-to-monopole ratio (although Fig. 2 shows that it does match the full 
^(a, tt) data very well). 



The presence of bias is an inevitable consequence of the nonlinear nature of galaxy for- 
mation, and the relation between mass and galaxy tracers is complex-'^^'-'^^''^^. However, 
there are good theoretical reasons to expect that h can indeed be treated as a constant 
on large scales, where the density fiuctuations are linear'^-'^''^'^. Redshift-space distortions 
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have thus been seen as an important method for weighing the universe'^^'^'*. To date, a 
number of papers have made significant detections of the Kaiser effect but the 
2dFGRS is the first survey that is large enough for the effect to be studied in detail. 

In order to estimate ^(u, tt), we follow standard methods'^'^''^^ that compare the 
observed count of galaxy pairs with the count estimated using a random distribution 
that obeys the same selection effects in redshift and sky position. These selection effects 
are well defined, but complex: the survey is tessellated into a pattern of 'sectors' defined 
by the overlap of the 2° diameter survey tiles, whose positions are chosen adaptively with 
the aim of being able to place a fibre on > 95% of the galaxies in the input catalogue. 
At the present intermediate stage of the survey, many tiles remain to be observed, and 
some regions of the survey presently contain redshifts for < 50% of the galaxies. Fur- 
thermore, the spectroscopic success rate (redshifts per allocated fibre) is > 95% in good 
conditions, but can fall to ~ 80% in marginal weather. We have implemented a number 
of independent algorithms for estimating the resulting survey selection effects, and are 
confident that we can measure the galaxy correlations robustly out to a separation of 
25h~^M.pc. For example, the redshift distribution in sectors with low spectroscopic 
completeness is biased to low redshifts, but it makes no significant difference whether 
or not we correct for this, or indeed whether the low- completeness regions are simply 
excised. In addition to allowing for survey completeness, it is necessary to give higher 
weight to regions with a low sampling density, to achieve the optimum balance between 
cosmic variance and shot noise^. In practice, we have chosen to truncate the analysis 
at a maximum redshift of ^ = 0.25. Within this volume, the exact optimum weight per 
galaxy varies very nearly as the reciprocal of the number density, so that all volume 
elements receive approximately equal weight. The redshift-space correlation function 
for the 2dFGRS computed in this way is shown in Fig. 2. The correlation-function 
results display very clearly the two signatures of redshift-space distortions discussed 
above. The 'fingers of God' from small-scale random velocities are very clear, as indeed 
has been the case from the first redshift surveys^. However, this is the first time that 
the detailed signature of large-scale flattening from coherent infall has been seen with 
high signal-to-noise. 

Quantifying redshift-space distortions 

The large-scale flattening of the correlation function may be quantifled by measuring 
the quadrupole moment of ^(cj, tt) as a function of radius. A negative quadrupole mo- 
ment implies flattening, whereas the flnger-of-God distortion tends to yield a positive 
quadrupole moment. Fig. 3 shows that the quadrupole-to-monopole ratio is positive 
on small scales, but that it falls with separation, becoming progressively more negative 
up to the largest separations at which it can be reliably measured. This arises partly 
because the underlying power spectrum is not a simple power law function of scale, so 
that the peculiar velocities have a different effect at different radii. By integrating over 
the correlation function, it is possible to construct quantities in which this effect is elim- 
inated. We shall not do this here, firstly because it seems desirable to keep the initial 
analysis as direct as possible. More importantly, finger-of-God smearing is a significant 
correction that will also cause the fiattening to depend on radius. We therefore have 
to fit the data with a two-parameter model, described in the caption to Fig. 2. The 
parameters are /3 and a measure of the size of the random dispersion in the relative 
velocities of galaxies, ap. In practice, ap plays the role of an empirical fitting parameter 
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to describe the scale on which the distortions approach the hnear-theory predictions. It 
therefore also incorporates other possible effects, such as a scale dependence of bias. 




r / h ^Mpc 



Figure 3 The flattening of the redshift-space correlation function is quanti- 
fied by the quadrupolc-to-monopole ratio, ^2/^0- This quantity is positive where 
fingers-of-God distortion dominates, and is negative where coherent infall dom- 
inates. The solid lines show model predictions for /? = 0.3 0.4 and 0.5, with a 
pairwise velocity dispersion of ap = 400kms~^ (solid lines), plus /? = 0.4 with 
dp = 300 and 500kms~^ (dashed lines). The ^2/^0 ratio becomes more negative 
as (3 increases and as Up decreases. At large radii, the effects of fingers-of-God 
become relatively small, and values of /? ~ 0.4 are clearly appropriate. 

The multipole moments of the correlation function are defined as ^e{r) = 
{2£ + l)/2 J_^^{(J — rsiiiO, IT — rcos9) P^(cos6') dcos9. In linear theory, the 
quadrupole-to-monopole ratio is given^^ by ^2/^0 = fi^) (4/3/3 + 4/3^/7) / (1 + 
2/3/3 -|- /3^/5). Here /(n) = (3 -|- n)/n, where n is the power-spectrum index of 
the density fluctuations: ^ oc r~^-^^'^K In practice, nonlinear effects mean that 
this ratio is a function of scale. We model this by using the real-space correla- 
tion function estimated from the APM survey-*^^'-*^®, plus the model for nonlinear 
finger-of-God smearing given in the caption to Fig. 2. 



The results for the quadrupole-to-monopolc ratio are shown in Fig. 3, which 
shows the average of the estimates for the NGP and SGP slices. The difference between 
the NGP and SGP allows an estimate of the errors to be made: these slices are inde- 
pendent samples for the present analysis of clustering on relatively small scales. For 
model fitting, it is necessary to know the correlation between the values at different r. A 
simple way of addressing this is to determine the effective number of degrees of freedom 
from the value of for the best-fitting model. A more sophisticated approach is to 
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Figure 4 Likelihood contours for f3 and the fingers-of-God smearing param- 
eter cTp, based on the data in Fig. 3 (considering 8h~^ Mpc < r < 25 Mpc). 
These are plotted at the usual positions for one-parameter confidence of 68% 
(shaded region), and two-parameter confidence of 68%, 95% and 99% (i.e. Ax^ = 
1, 2.3, 6.0, 9.2). The maximum-likelihood solution is (3 — 0.43 and dp = 385kms~^ 
The value for the large-scale pairwise dispersion is in reasonable agreement with 
previously suggested values^^; however, for the present analysis CTp is an unin- 
teresting parameter. If we marginalize over dp (i.e. integrate over Up, treating 
the likelihood as a probability distribution), the final estimate of (3 and its rms 
uncertainty is /? = 0.43 ± 0.07. 

We believe that this result is robust, in the sense that systematic errors 
in the modelling are smaller than the random errors. We have tried assuming 
that the power spectrum for k < 0.1hMpc~^ has the shape of a = 0.3 ACDM 
model, rather than the APM measurement; this has a very small effect. A more 
serious issue is whether the pairwise velocity dispersion of galaxies may depend 
strongly on separation, as is found for mass particles in numerical simulations^ 
Assuming that the pairwise velocity dispersion cTp rises to twice its large-scale 
value below 1 Mpc reduces the best-fit f3 by 0.04. This correction is small 
because our analysis excludes the nonlinear data at r < 8 Mpc. 



generate realizations of ^ (cr, tt) , and construct the required covariance matrix directly. 
One way of achieving this is to analyze large numbers of mock surveys drawn from 
numerical simulations^"^. A more convenient method is to generate direct realizations 
of the redshift-space power spectrum, using Gaussian fiuctuations on large scales, but 
allowing for enhanced variance in power on non-linear scales^^'^^'^^. In practice, the 
likelihood contours resulting from this approach agree well with those from the simple 
approach, and we are confident that the resulting errors on /3 are realistic. These con- 
tours are shown in Fig. 4, and show that there is a degree of correlation between the 
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preferred values of (3 and Up, as expected. For present purposes, Up is an uninteresting 
parameter, so we marginalize over it to obtain the following estimate of j3 and its rms 
uncertainty: 

/? = 0.43 ± 0.07. 

This result is the first precise determination of fi from redshift-space distortions. 
The best previous studies-*^®'^^'^® have achieved no more than a detection of the effect 
at the 3(7 level. Before discussing the implications of our result, we should therefore 
consider some possible small systematic corrections that have been unimportant in 
earlier work. First, the Kaiser analysis applies only in the small-angle approximation, 
and in principle corrections might be needed for wide-angle surveys such as ours^^. 
However, with our weighting scheme, the mean angular separation of pairs with spatial 
separations < 30 /i"-*^ Mpc is only 2.5°, so this is not a concern. There is potentially 
a significant correction for luminosity effects. The optimal weighting means that our 
mean luminosity is high: it is approximately M^j = —20.3, or 1.9 times the characteristic 
luminosity, L*, of the overall galaxy population^'^. A number of studies^^ have suggested 
that the strength of galaxy clustering increases with luminosity, with an effective bias 
that can be fitted by h/h* = 0.7 + 0.3(L/L*). This effect has been controversial^^, 
but the 2dFGRS dataset favours a similar luminosity dependence, as will be described 
elsewhere. We therefore expect that (3 for L* galaxies will exceed our directly measured 
figure. Applying a correction using the given formula for b{L), we deduce 

f3(L = L*) = 0.54 ±0.09. 

Finally, the 2dFGRS has a median redshift of 0.11. With weighting, the mean redshift 
in the present analysis is ^ = 0.17, and our measurement should be interpreted as j3 
at this epoch. The extrapolation to 2; = is model-dependent, but probably does not 
introduce a significant change^°. 

Consistency with microwave-background anisotropies 

Our results are significant in a number of ways. First, we have verified in some detail 
that the pattern of redshift-space distortions associated with the gravitationally-driven 
growth of clustering exists as predicted. Although gravitational instability is well estab- 
lished as the standard model for the formation of large-scale structure, it is an important 
landmark to have verified such a characteristic feature of the theory. Extracting the full 
cosmological implications of our measurement of requires us to know the bias pa- 

rameter in order to determine O. For example, our measurement implies Q, = 0.36 ±0.10 
if L* galaxies are imbiased, but it is difficult to justify such an assumption. In principle, 
the details of the clustering pattern in the nonlinear regime allow the — 6 degeneracy 
to be broken, yielding a direct determination of the degree of bias^^. We expect to 
pursue this approach using the 2dFGRS. For the present, however, it is interesting to 
use an independent approach. Observations of anisotropies in the cosmic microwave 
background (CMB) can in principle measure almost all the cosmological parameters, 
and current small-scale anisotropy results are starting to tighten the constraints. In a 
recent analysis"^"*", best-fitting values for the densities in coUisionless matter (c), baryons 
(6), and vacuum (v) have been obtained: Qc + ^b + ^v = 1-11 ±0.07, Qch^ = 0.14 ±0.06, 
ribh"^ = 0.032±0.005, together with a power-spectrum index n = 1.01±0.09. Our result 
for P gives an independent test of this picture, as follows. 



8 




0.01 0.1 1 

k / h Mpc"^ 



Figure 5 The dimensionlcss matter power spectrum at zero redshift, A^(/c), 
as predicted from the aUowed range of models that fit the microwave-background 
anisotropy data, plus the assumption that Hq = 70 kms^^Mpc"-*^ ± 10%. The 
solid line shows the best-fit model^-*- (power-spectrum index n = 1.01, and density 
parameters in baryons, CDM, and vacuum of respectively 0.065, 0.285, 0.760). 
The effects of nonlinear evolution have been included'*^ . The shaded band shows 
the la variation around this model allowed by the CMB data. The solid points 
are the real-space power spectrum measured for APM galaxies. The clear conclu- 
sion is that APM galaxies are consistent with being essentially unbiased tracers 
of the mass on large scales. Since the CMB data also constrain the range of 
this allows P to be predicted. 



The only parameter left undetermined by the CMB data is the Hubble constant, 
h. Recent work^^'^^ indicates that this is now determined to an rms accuracy of 10%, 
and we adopt a central value oi h = 0.70. This completes the cosmological model, 
requiring a total matter density parameter Q = Qc + = 0.35 ± 0.14. It is then 
possible to use the parameter limits from the CMB to predict a conservative range for 
the mass power spectrum at z = 0, which is shown in Fig. 5. A remarkable feature of 
this plot is that the mass power spectrum appears to be in good agreement with the 
clustering observed in the APM survey. For each model allowed by the CMB, we can 
predict both b (from the ratio of galaxy and mass spectra) and also P (since a given CMB 
model specifies O) . In practice, we determine b by determining the mean ratio of power 
spectra over the range 0.02 < k < 0.1 /iMpc~^, where the APM measurement is robust 
and where scale-dependent bias and nonlinearities should be unimportant. Considering 
the allowed range of models, we then obtain the prediction /3cmb+apm = 0.57 ± 0.17. A 
flux- limited survey such as the APM will have a mean luminosity close to L*, so the 
appropriate comparison is with the 2dFGRS corrected figure of = 0.54 ± 0.09 for L* 
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galaxies. These numbers are in very close agreement. 



This analysis of galaxy clustering in the 2dFGRS thus gives strong support to the 
simplest picture of cosmological structure formation, in which the primary mechanism 
is gravitational instability in a sea of coUisionless dark matter. We have shown that the 
fluctuations seen in the CMB (which measure structure at a redshift z 1100) can be 
extrapolated to the present to predict the peculiar velocities that distort redshift-space 
clustering. The agreement between this extrapolation and direct observations from the 
2dFGRS is a remarkable and highly non-trivial test of the basic model. The precision of 
data in both areas should improve rapidly, and the use of /3 as a meeting ground between 
studies of the CMB and large-scale structure will undoubtedly lead to more demanding 
tests of the theory in years to come. For the present, we can say that there is complete 
consistency between clustering in the 2dFGRS and the emerging 'standard model' of 
cosmology: a spatially flat, vacuum-dominated universe with density parameter Q ~ 0.3. 
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